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Type-I proton superconducting cores of neutron stars break up in a magnetic field into alternating 
domains of superconducting and normal fluids. We examine two channels of superfluid-normal fluid 
friction where (i) rotational vortices are decoupled from the non-superconducting domains and the 
interaction is due to the strong force between protons and neutrons; (ii) the non-superconducting 
domains are dynamically coupled to the vortices and the vortex motion generates transverse electric 
fields within them, causing electronic current flow and Ohmic dissipation. The obtained dissipation 
coefficients are consistent with the Eulerian precession of neutron stars. 
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I. INTRODUCTION 



The timing observations of radio-pulsars provide a 
unique tool to study the properties of superdense matter 
in compact stars. While pulsars are known to be prefect 
clocks over long periods of time, the timing observations 
of past few decades revealed several types of timing "ir- 
regularities" in a subclass of isolated compact objects. 
Glitches (or macrojumps) - sudden increases in the pul- 
sar rotation frequency and its derivative - are the most 
spectacular examples of timing anomalies. The slow re- 
laxation of their spin and its derivative following a glitch 
has been interpreted as an evidence for superfluidity of 
compact star interiors [l| . At temperatures prevailing in 
an evolved compact star the dense hadronic matter is ex- 
pected to be in the superfluid state due to the attractive 
component of the nuclear force which binds neutrons into 
Cooper pairs either in the relative 1 5'o state (at low den- 
sities) or 3 P 2 — 3 F 2 state (at high densities) 0. Other 
members of the spin-1/2 octet of baryons (protons and 
strangeness 5=1 particles £ _0+ , A) are expected to 
pair in the 5-wave channel due to their low concentra- 
tion. If the central densities of compact stars exceed 
the density of deconfincmcnt phase transition to a quark 
matter phase, the deconfined quark matter will be in one 
of the many possible color superconducting states Q. 
Since the glitches are related to the axisymmetric per- 
turbations from the state of uniform rotation there is a 
twofold degeneracy in the interpretation of the data on 
both the jumps and the post-jump relaxations; the time- 
scales of these processes can be associated either with 
the weak or strong coupling between the superfluid and 
normal fluid in the star's interiors, and it is impossible to 
distinguish between these regimes on the basis of glitch 
observations alone. During the recent years it became 
increasingly clear that another type of timing anomaly 
- the long term periodic variations superimposed on the 
spin-down of the star - can provide an additional and 
independent information on the dynamical coupling be- 
tween the superfluid and the normal fluid in compact 
stars. If these irregularities are interpreted in terms of 
the precession of the star (a motion which involves non- 
axisymmetric perturbations from the rotational state) 
the degeneracy inherent to the interpretation of glitches 



is lifted. It turns out that the free precession is pos- 
sible only in the weak coupling limit and it is damped 
in the strong coupling case. Interpretations of the tim- 
ing anomalies in pulsars in terms of the friction between 
the superfluid and the normal fluid require a model of 
the friction between the superfluid and normal compo- 
nents of the star on mesoscopic scales characteristic for 
the vorticity. This paper discusses two new mechanism 
of mutual friction in the core of a neutron star in the 
case where protons form a type-I superconductor. The 
remainder of the introduction sets the stage by briefly 
reviewing the relevant physics. Section [Q] studies the dy- 
namics of a type-I superconducting model where there is 
a single normal domain per rotational vortex. In section 
IIIII we discuss the dissipation in an alternative picture 
where there is a large number of rotational vortices asso- 
ciated with a single normal domain. Section llVl is devoted 
to the implications of the dissipative dynamics of type-I 
superconductors for the free precession of compact stars 
and contains a brief summary of the results. 



A. No-go theorems for precession 

To see how the superfluidity of neutron stars changes 
their Eulerian precession (which would be intact if the 
neutron stars were non-superfluid) let us begin with the 
equation of motion of approximately massless neutron 
vortex 



Psk (vs — vl) x n — rju — rj ' [n X tt) = 0, 



(1) 



where u = % - «l, and v$, vjy and vl are the veloc- 
ities of the superfluid, the normal fluid and the vortex; 
ps is the effective neutron density, k is the quantum of 
circulation, n = k/k, and the coefficients r\ and rf are the 
measure of the friction between the neutron vortex and 
the ambient normal fluid. Here we work within the two- 
fluid superfluid hydrodynamics, where it is assumed that 
the hydrodynamic forces are linear functions of the veloc- 
ities, which guarantees that the energy variation is always 
a quadratic form. Shaham first observed that the long- 
term, Eulerian precession is impossible if the neutron 
vortices are strongly pinned In terms of the meso- 
scopic parameters in Eq. QJ, his observation is equiv- 
alent to the statement that precession is absent in the 
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limit £ — > co, £' - ► where £ = tj/psk and £' = rj'/psK 
are the drag-to-lift ratios. Since in the frictionless limit 
a star must precess at the classical frequency eSl, where 
e is the eccentricity and Q is the rotation frequency, it 
is clear that there exists a crossover from the damped to 
the free precession as C is decreased. The crossover is de- 
termined by the dimensionless parameters (Is /In) and 
(Is /In) 0' where Is is the moment of inertia of the su- 
perfluid and In is the moment of inertia of the crust plus 
any component coupled to it on time-scales much shorter 
than the precession timescale and = (/[(I - C') 2 + C 2 
j3' = 1 — 0(1 — C')/C- The precession frequency is 
(hereafter SWC) 



1 1 

5] 



flp = efls 



1 + 0' 



i0- 



(2) 



where Qs is the spin frequency and e is the eccentric- 
ity. The result of SWC can be cast in a no-go theo- 
rem that states that the Eulerian precession in a super- 
fluid neutron star is impossible if (Is /In) C > 1 (assum- 
ing as before £' — > 0). There is a subtlety to this re- 
sult: the precession is impossible because the precession 
mode, apart from being damped, is renormalized by the 
non-dissipative component of superfluid-normal fluid in- 
teraction (oc 0'). In effect the value of the precession 
eigen-frequency drops below the damping frequency for 
any £ larger than the crossover value. Note that this 
counter-intuitive result can not be obtained from the ar- 
guments based solely on dissipation: in fact, according to 
Eq. (J2J the damping time-scale for precession increases 
linearly with £ and in the limit ( — *■ oo one would pre- 
dict (wrongly) undamped precession. If a neutron star 
contains multiple layers of superfluids the picture is more 
complex, but the generic features of the crossover are the 
same 



B. Previous work 

Long term variabilities were observed in a number of 
pulsars and have been attributed phenomenologically to 
precession of the neutron star (see ref. 6] and refer- 
ences therein). A strong case for long-term variability 
(again attributed to precession) was made recently by 
Stairs et al Q. While it is common to study pertur- 
bations from the state of uniform rotation, Wasserman 
[U demonstrated that the precessional state may corre- 
spond to the local energy minimum of an inclined rotator 
if there is a large enough magnetic stress on the star's 
core. This type of precession is likely to be damped away 
by the superfluid-normal friction. Link [j| argued that 
the long-term variations, which can be fitted by assum- 
ing Eulerian precession of the pulsar, are incompatible 
with type-II superconductivity of neutron stars. Type-I 
superconductivity was proposed to resolve the discrep- 
ancy 0. Jones ^3] argued that the friction of vortices 
in the crusts against the nuclear lattice will give rise to 
a dissipation which will damp the free precession; thus, 



the free precession (even in absence of pinning) would 
be incompatible with the known properties of matter at 
subnuclear densities. 



C. Type-I superconducting neutron stars 

As is well known, type-I superconductivity arises when 
the Ginzburg-Landau parameter satisfies the condition 
kgl = V£ — l/v2) where A is the magnetic field pene- 
tration depth and £ is the coherence length. Type-I su- 
perconductivity can arise loc ally within the current mod- 
els based on the BCS theory [Tj, with domain structures 
analogous to those observed in laboratory experiments. 
In ref. [Tl| the theory of these structures was constructed 
along the lines of the theories developed for laboratory 
superconductors, where the magnetic fields are generated 
by normal currents driven around a cylindrical cavity by 
temperature gradients However, global type-I su- 

perconductivity would require a suppression of the pro- 
ton pairing gap A p (due to the scaling kql oc £ _1 oc A p ) 
by polarization or related effects. Buckley et al stud- 
ied the effect the interactions between the neutron and 
proton Cooper pairs would have on the type of the pro- 
ton superconductivity. Their results suggest that type-I 
superconductivity can be enforced within the entire core 
without the suppression of the pairing gap if the strength 
of the yet unknown interaction between Cooper pairs will 
turn out to be significant. 

The equilibrium structure of the alternating supercon- 
ducting and normal domains in a type-I superconduc- 
tor is a complicated problem and depends, among other 
things, on the nucleation history of the superfluid phase. 
The equilibrium dimension of a layer is of the order of 
magnitude d ~ yL^ where L is the size of the core; (for 
typical parameter values L ~ 5 x 10 5 cm and £ ~ 200 fm, 
d ~ 3.2 x 10~ 3 cm). By flux conservation, the ratio of the 
sizes of the superfluid and normal domains is given by the 
relation ds/dN = y/H cm /B ~ 10, where B ~ 10 12 G is 
the average value of the magnetic induction, H cm ~ 10 14 
G is the thermodynamic magnetic field. The dimensional 
analysis above suggests that there is roughly a single 
normal domain per neutron vortex. We shall consider 
below the dynamics of two distinct models where (i) a 
neutron vortex features a single coaxial normal domain 
of a smaller size according to ref. [Tll| (hereafter 1 + 1 
model) and (ii) a large fraction of neutron vortices is ac- 
commodated by a single normal domain, as described in 
ref. [lj| (1 + N model). The difference between these 
models is not simply the number of the neutron vortices 
accommodated by a single non-superconducting domain; 
in the 1+1 model the magnetic fields are generated dy- 
namically (this is explained in more detail below) and 
hence the normal domains are tied to the neutron vor- 
tices on dynamical times scales. In effect the motion of a 
neutron vortex requires the motion of the normal domain 
attached to it on the dynamical timescale. Therefore the 
dissipation is due to the interaction between the com- 



3 



bined structure (a vortex plus a normal domain) with 
the background electron liquid. While it is conceivable 
that the normal domains are tied up to the macroscopic 
neutron vortex lattice on certain timescales in the 1 + N 
model, we shall assume that the neutron vortex lattice 
sweeps through the non-superconducting domains and 
the dissipation arises from the scattering of the normal 
protons off the cores of rotational vortices 01 ■ 




II. DYNAMICS OF 1+1 STRUCTURES 

In two component superfluids the supercurrent of any 
given component transports mass of both components; 
this is the essence of the cntrainmcnt effect first studied in 
the context of charge neutral and non-rotating superfluid 
mixtures of 3 He- 4 He ^BJ. If one of the components is 
charged - a case first studied in refs. 0,EJ - the neutral 
superfluid which rotates by forming a lattice of charge- 
neutral vortices generates magnetic fields because neutral 
supercurrent carries along a finite mass of the charged 
component. The mass currents of neutrons and protons 
Pp/n (here and below indices p and n refer to the protons 
and neutrons), are related to their velocities v„, and v p 
by a density matrix in the isospin space [TfJ H3 



Pn 



Ppp Ppn 
Pnp Pnn 



(3) 



where in the mean-field approximation the elements of 
the density matrix can be expressed through the effec- 
tive masses of neutrons and protons. Mendell obtained 
previously the general form of mutual friction damping 
from vortices, which incorporated the entrainment effect, 
in the case of type-II superconductivity |l8| . 

If the proton superconductor is type-I, non- 
superconducting domains coaxial with the neutron 
vortices nucleate in response to the entrainment current 
set-up by the vortex circulation Consider a cylin- 

drical domain of radius a coaxial with a vortex (see Fig. 
1 for an illustration). In the cylindrical coordinates (r, 
<f>, z) with the symmetry axis at the center of the vortex 
the magnetic field induction is [ref. [Tl|, Eq. (20)] 



$1 Tn ( b - 



2ttA 2 



N(r) 



a ) N(a) ' 



B X = By — 0, 



(4) 



where $i = (p pn / p PP )$o , &o is the flux quantum, A 
is the magnetic field penetration depth, b is the vortex 
(outer) radius and 



JV(r) = Jo 



K 



(5) 



where Iq(z) and Kq(z) are the modified Bessel functions. 
In the mean-field approximation the magnitude of the 
(non-quantized) flux $i is determined by the effective 
mass of a proton quasiparticle p pn / ' p pp = \m*/m p — 1|. 

Now we can write down the proton supercurrent j<f, = 
(c/47r)(V x B)^ by substituting the S-field from Eq. 




FIG. 1: An illustration of the structure of a rotational vortex 
placed in a type-I superconductor. The vortex velocity field 
is shown by the concentric circles. The non-superconducting 
domain (shaded region) of radius a is coaxial with the vortex 
and carries a magnetic field H C m ~ 10 14 G. The vortex motion 
along the :r-axis generates a transverse electric field, which 
drives the electron current through the domain and leads to 
Ohmic dissipation. 



Q. We need, however, the velocity of superfluid pro- 
tons, which is the difference between the net supercurrent 
and the supercurrent moving with the neutron superfluid 
velocity; we find 



kh 
2mr 







[H 









coth 



f 



(6) 



Eq. © keeps the leading order term of the expansion of 
the Bessel functions with respect to large arguments r/X, 
where r is a mesoscopic scale > a. 

Consider now a vortex which moves at a constant ve- 
locity vl, and carries a coaxial normal domain of protonic 
fluid with respect to the background electron liquid (see 
Fig. 1). The equation of motion of the superfluid pro- 
tons, written in the reference frame where the vortex is at 
rest, acquires an additional term (vl ■ V) • v p . The conti- 
nuity of the electro-chemical potentials of the superfluid 
and normal phases across the boundary of the normal 
and superconducting phases gives ps — Pn, where 



y 

p s = p- mlv L v p (r)-, 
1 r 



Pn = fi-eq 



(7) 



here p is the proton chemical potential in equilibrium, 
4> is the scalar electric potential and we use 2d Carte- 
sian coordinates with vortex circulation along z-axis and 
vortex velocity along the x-axis, see Fig. 1; small terms 
0(A 2 /p) are neglected. Thus, the motion of the vortex 
generates a constant transverse electric field across the 
normal domain 



Ey = -\7y(j) = 



m*v L v p (a) 



E T — E, 



0. 



(8) 



The power dissipated per unit length of a vortex is W = 
aE 2 (a/b) 2 , where a is the electrical conductivity and the 
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TABLE I: Listed are the protons density (column 1), the 
Fermi- wave number (2) , the effective mass of protons (3) , the 
magnetic field penetration depth (4), the critical thermody- 
namic field (5), the electron relaxation time (6), the electrical 
conductivity (7), and the drag-to- lift ratio (8). The short- 
hand notions pi4 = p/10 14 , etc., are used. 
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factor (a/b) 2 is the fractional 2d volume occupied of the 
domain. Upon substituting Eq. (jSJl in this relation, we 
obtain an alternative form of dissipation W — r\v\ , which 
identifies the friction coefficient 



1 c 2 \2irabJ 



a fb\ , (b 
—In — ] coth 
A V a 



Equation @ is our central result, which defines the fric- 
tion coefficient for a vortex featuring a single coaxial do- 
main of a type-I proton superconductor in terms of the 
electrical conductivity of the electron Fermi- liquid in nor- 
mal proton matter. 

Since the mean-free path of electrons is much smaller 
than the size of a single normal domain, wc can neglect 
the finite-size effects and use the result for the bulk nor- 
mal matter |19). The zero-field conductivity of ultra- 
relativistic electrons is (J = n e e 2 CT c f '(hkp), where the 
relaxation time for the Coulomb scattering of electrons 
off the protons in the normal domains is LL9| 



12 

72 



kFT 

ckp 



(10) 



where ep and kp are the electron Fermi-energy and 
Fermi-wavenumber, and kpp = [Akpmpe 2 /^ 2 ] 1 / 2 is the 
Thomas-Fermi wave-length. Since the Larmor radius of 
an electron moving in a normal domain carrying a field 
ff cm ~ 10 14 G is much smaller than the linear size of the 
domain, the conductivity becomes a = cto/(w c t c ) 2 , where 
lj c = (eH cm ) I '(hkp), where lj c is the electron cyclotron 
frequency. Table I lists some of the relevant parameters 
of the problem, including transport quantities and the 
drag-to-lift ratio at the temperature T = 10 s K. The ki- 
netic coefficients can be rescaled to other temperatures 
by using the scalings r oc T~ 2 , a oc T 2 and £ oc T 2 . The 
drag-to-lift ratio satisfies the condition £ <C 1 at all tem- 
peratures (the largest values correspond to temperatures 
just below the critical temperature T c ~ 10 9 K). We shall 
return to the implications of these results for the neutron 
star precession in the closing section. 



III. DYNAMICS OF 1+N STRUCTURES 

Buckley et al |l| argued (qualitatively) that the size of 
the normal domains could be large enough to accommo- 
date about N = 10 neutron vortices across a single nor- 
mal domain of protonic fluid. Since there is no dynam- 
ical coupling (in the sense of the entrainment) between 
the vortices and the normal domains the damping of the 
differential rotation between electron-proton plasma and 
the neutron supcrfluid is due to the interaction of domain 
(non-superconducting) protons with the core quasiparti- 
cles confined in the neutron vortex core. The relaxation 
process is thus the same as for the case where the proton 
fluid is non-superconducting over the entire bulk of the 
core, but the final result needs to be rescaled by the ratio 
of the areas occupied by the normal and superconducting 
layers. The relaxation time per single vortex is |14| 



G 



kFn 



m n fi 



pn 



hm*Ta np 



exp 0.02 



£FnT 



(11) 



(9) where kp p and kp n are the Fermi-wave-numbers of pro- 
tons and neutrons, /i* n = TO*m*/(m* + m*) is the re- 
duced effective mass, with m* begin the neutron effective 
mass, a is the total in-medium neutron-proton scattering 
cross-section, A n is the gap in the neutron quasiparticle 
spectrum, ep n is the neutron Fermi-energy. [Eq. i|ll|) 
differs from the analogous expression in ref. |l4| by the 
factor Am n /hP; here P is the pulsar period, m n - free- 
space neutron mass]. 

In the relaxation time-approximation, the force exerted 
by the normal proton on a single vortex is given by a 
phase-space integral 



/ = — 

Trvn 



dp 



/f(p,«) = T)U, 



(12) 



over the proton Fermi-distribution function, fp(p,u) = 
[exp(e p — ep p + p ■ u)/T + 1] _1 , where the quasiparticle 
energy shifted due to the motion with a velocity u. Here 
e p is the dispersion relation of normal protons, ep p is 
their Fermi-energy. The integral l|12|l is straightforward 
in the T = limit and we obtain 

hkppUp 



V 



(13) 



1 up 



where n p is the proton number density. The result for 
the friction coefficient r\ and the corresponding drag-to- 
lift ratio for several densities are listed in Table II for the 
case where the proton fluid is non-superconducting. For 
a given model of the type-I superconducting structure, 
the friction coefficient r\ must be rescaled by a factor 
(d N /d s ) 2 . 



IV. IMPLICATIONS FOR PRECESSION AND 
CONCLUSIONS 

The SWC no-go theorem requires the condition 
(Is/In)( < 1 to be fulfilled for precession to occur; 
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TABLE II: The columns 3-5 list the relaxation time for proton 
scattering off the neutron quasiparticles in the neutron vortex 
cores (3), the vortex friction coefficient (4), and the drag-to- 
lift ratio (5) for the temperature T — 10 8 K. The columns 7-9 
list the same parameters for T = 10 7 K. 
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otherwise the precession is damped. The magnitude of 
the ratio Is /In depends on the superfluid-normal fluid 
friction within all superfluid regions of a neutron star 
and is difficult to access. Glitches and post-glitch re- 
laxation provide a model independent lower bound on 
Is /In > 0.1. An upper bound is difficult to place, since 
the deep interiors of neutrons stars, if superfluid, could 
be decoupled from the observable parts of the star on 
evolutionary timescales without any effect on short time- 
scale physics |23 (but one needs ( — > 0, rather than 
£ — > 00, to prevent the damping of the precession). How- 
ever, it is rather unlikely that this ratio exceeds unity by 
many orders of magnitude. For the first dissipation chan- 
nel studied C ~ 1CT 4 - 1(T 3 (Table I) and this clearly 
suggests an undamped precession. The second chanel is 
more effective, £ ~ 1CP 2 (Table II), but these numbers 
must be reduce by a factor (ds /djv) 2 — 100 . On account 
of the lower bound on the ratio of the moments of iner- 
tia, one can conclude that the precession is undamped 
for both dissipation mechanism. 

We have provided a first discussion of the dynamics 
of the type-I superconducting domains in neutron star 
interiors. Although the details of the coupling of the 



superconducting and normal components of the star de- 
pend on the form of the superconducting-normal struc- 
tures that nucleate and, in particular, whether the non- 
superconducting domains are dynamically coupled to the 
rotational vortices or not, in all cases we find friction co- 
efficients that imply undamped precession. 

The results above by no means suggest that type-I su- 
perconductivity is the only resolution to the precession 
puzzle and alternatives should be searched for. An alter- 
native to free, Eulerian precession is the forced precession 
due to time-dependent periodic torques Q. Other peri- 
odic motions, for example, Tkachenko oscillations of the 
vortex lattice could generate the observed timing features 

mmm. whnethe 

eigen- frequencies of the Tkachenko 
modes are of correct order of magnitude and could ex- 
plain long-term periodicities it remains to be studied 
whether these modes will be undamped by the mutual 
friction between the superfluid and the normal fluid. 

The fact that statistically insignificant number of pul- 
sars show long-term variabilities, indicates that a subtle 
tuning is needed for the underlying mechanisms to work. 
On the other hand, the Eulerian precession, if undamped, 
should be a common place in the pulsar population, since 
neutron stars frequently undergo non-axisymmetric per- 
turbations such as glitches and quakes. 
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